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Abstract 

Particle sedimentation in the vicinity of a fixed horizontal vortex with time-dependent intensity 
can be chaotic, provided gravity is sufficient to displace the particle cloud whilst the vortex is 
off or weak. This "stretch, sediment &; fold" mechanism is close to the so-called blinking vortex 
effect, which is responsible for chaotic transport of perfect tracers, except that in the present case 
the vortex motion is replaced by gravitational settling. In the present work this phenomenon is 
analyzed for heavy Stokes particles moving under the sole effect of gravity and of a linear drag. 
The vortex is taken to be a fixed isolated point vortex the intensity of which varies under the effect 
of either boundary conditions or volume force. When the unsteadiness of the vortex is weak and 
the free-fall velocity is of the order of the fluid velocity, and the particle response time is small, 
the particle motion equation can be written asymptotically as a perturbed hamiltonian system the 
phase portrait of which displays a homoclinic trajectory. A homoclinic bifurcation is therefore likely 
to occur, and the contribution of particle inertia to the occurrence of this bifurcation is analyzed 
asymptotically by using Melnikov's method. 

Key-words : particle-laden flows, inertial particles, sedimentation, chaotic motion, homoclinic 
bifurcation. 

1 Introduction 

In contrast with the chaotic advection of perfect tracers, which has been a topic of great interest 
in the last decades (Arnold j^, Henon [9], Aref |3], Ottino |13)). little is known about the chaotic 
motion of inertial particles in laminar flows. Chaotic advection is entirely contained within an 
elementary equation of kinematics, that is : 

where ^(t) is the position of a perfect tracer advected by the velocity field f- This velocity field 
is not required to be turbulent for Lagrangian chaos to take place : it can be laminar and entirely 
specified. In particular, chaotic advection provides efficient mixing properties, and this is one of the 
reasons why this effect attracted the attention of the physics community. Mixing is indeed a key 
phenomenon in nature, and a challenge for engineers as soon as one has to mix large amounts of very 
viscous fluids, or small quantities of fluid in tiny domains. Nevertheless, the problem "starts rather 
than ends with the speciflcation of the velocity field" (Ottino |14|). so that many analyses are still 
being performed to quantify the Lagrangian properties of various flows. Among the mechanisms 
leading to chaotic advection, the blinking vortex is probably the easiest to perform (Aref [2]). Two 
key ingredients are required here : the differential rotation and the displacement of the vortex. The 
former ingredient is responsible for stretching and the latter is required for folding to take place. 
Many chaotic flows are based on this property and are described in several reviews [3]' |13) . 

In the case where the particle is not a pure tracer, surprising chaotic trajectories can also occur, 
even when the inclusion obeys a linear drag law. For example, aerosols can have complex trajectories 
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in elementary cellular or ABC flows (Maxey & Corrsin |12) ; Wang , Maxey, Burton & Stock |19| ; 
Mac Laughlin |10| ; Fung [7] ; Tsega, Michaelides &; Eschenazi [T7] ; Rubin, Jones & Maxey |16|). 
The complex trajectories reported in these works are due to the non-uniformity of the flow, to the 
finite response time of the inclusion, and to gravity if any. 

The goal of the present work is to show that heavy particles with small response time but 
non-negligible terminal velocities can have chaotic trajectories under the combined effect of a time- 
periodic differential rotation (which creates stretching of a particle cloud) and of gravity (which 
prepares the particle cloud to undergo folding). In particular, we will consider flows where un- 
steadiness does not suffice to induce Lagrangian chaos, like two-dimensional fiows of the form : 

]^f(jt,t) = \^^f(]t){l + esmut) (2) 

with '^^(jI) corresponding to a horizontal vortex. Indeed, one can check that the dynamical system 
([1]) is not chaotic in this case, and that fluid points go to and fro along the streamlines of the vortex. 
Clearly, unsteadiness does not produce chaos because the vortex does not move : not only should 
we switch the vortex off sometimes, but also should we light it up somewhere else to produce chaos 
according to the blinking vortex mechanism. However, if gravity is sufficient to displace the particle 
cloud whilst the vortex is off (or weak), one could expect chaos according to some kind of "gravity- 
induced blinking vortex". Note, however, that the key role of gravity in this scenario could also be 
played by electrostatic forces, or swimming (if particles are bacteria or plankton), to name but a 
few examples. 

In the following we show that heavy Stokes particles can undergo such a mechanism. The dynamics 
of inertial particles is much more complicated than the one of tracers, since one has to solve for 
both the equations of the flow induced by the inclusion, and the motion equations of the inclusion. 
Nevertheless, significant simplifications arise when the fiow induced by the inclusion is a quasi- 
steady creeping flow. In the case of tiny heavy particles (e.g. aerosols) carried by a fluid with 
infinite extent one often writes 

Xp = ^ + -(Vf{jlp,t)-Xp), (3) 

where ^p(t) is the particle position, Tp denotes its response time, and t is the terminal velocity of 
the inclusion in the very same fluid at rest. This is the simplest equation for non-ideal tracers, which 
requires the particle Stokes and Reynolds numbers to be much smaller than unity, and brownian 
diffusion to be negligible. In the present paper we will assume that Q is non-dimensionalized by 
the typical length scale and velocity of T^j respectively, and that : 

rp< 1 and Vt = | V^t| = 0(1)- 

The former condition can be thought of as a consequence of the fact that the viscous time scale 
over the particle radius is much smaller than convective flow time scales. The latter manifests 
non-negligible sedimentation effects. If in addition we assume e ^ 1, the particle motion equation 
contains two independent small parameters, namely e and r^. Classical asymptotic expansions 

of the form X„ -|- O(e^) enable one to write the particle motion equation as a 

three-dimensional non-autonomous dynamical system : 

?p = V^T + ^)(^p) + £(^}(^p, t) - k^^).(^) + V^t)) + O(e'), (4) 

where we have set Tp = ke, with k held flxed as e — )> 0. Because the flow is 2D and div (V^t + ^j) = 
0, Eq. d?]) is a perturbed hamiltonian system with one and a half degree-of-freedom. The phase 
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portrait of the unperturbed system (e = 0) can display homoclinic or heteroclinic trajectories, 
which are key-ingredients of chaos for such system. Figure [T]J^b) shows such trajectories for particles 
moving in the vicinity of a point vortex the streamlines of which are sketched in Fig. [1] (a). The 
homoclinic trajectory (dashed line) links the saddle point where the fluid velocity V^j(^p) balances 

the terminal velocity (Throughout this paper we have chosen = -VT~^y, Vt>0). Such 

particle trajectories have often been observed in particle-laden flows. For example, trajectories like 
Fig. mb) have been investigated by Davila & Hunt [6]. 

The perturbations contained in the 0(e) terms of Eq. (jj]) can have several tremendous effects on 
the particle dynamics. In particular, a homoclinic bifurcation can occur, leading to chaotic particle 
settling or trapping. The purpose of the present paper is to investigate under which conditions such 
a bifurcation could occur. 



2 Asymptotic analysis in the vicinity of the separatrix 



The basic vortical flow investigated in this paper is the point vortex : 

T^°.(^) = rn{r)'^g, with n{r) = Oo(r/iio)"^ 



(5) 



This velocity field is set non-dimensional in the following by using i?o as length units and l/r^o a-s 

time units. The streamfunction ^/^'^ of this non-dimensional flow therefore reads : ijj^(ip{t)) = ——Log\^p\'^. 

To leading order 0{e^) the particle dynamics reads : 

Xp=}^T + '^}(^p) (6) 

and the corresponding trajectories therefore correspond to the iso-values of the hamiltonian : 

H{x, y) = ip°{x, y) + xVt 

and are sketched in figure [1] (b). In order to investigate separatrix splitting under the effect of the 
O(e^) terms it is necessary to solve analytically the leading-order motion ([6]) with ^p(iboo) = A 
(saddle point). To our knowledge, even for the simple flow considered here, this cannot be done. 
To be precise one cannot obtain a simple solution which would make the Melnikov integral easy to 
calculate analytically. However, one can obtain semi-analytical results by rescaling the variables. 
Indeed, we set : 

^ (t) = 1-^{t) with r = tV^ 
Vt 

and notice that the velocity field investigated here satisfies : 



The particle dynamics therefore reads : 



d 



dT 



+ e 



(7) 



The leading order dynamical equation is now independent of Vt, and the corresponding phase 
portrait in the 1^ plane is similar to the one of figure [T](b), with the saddle point located at (1,0) 
(figure [SJa)). Let 1^o(t) be a solution of the leading-order dynamics : 



(1,0) 



(8) 
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(the dot upon a 1^ indicating a derivation with respect to r) which satisfies, accordingly : 



— > 



/•V " /v-f o; - e yj. 

When the 0(e) terms are taken into account, and because the perturbation of the rescaled system is 
time-periodic (with period Ty = TV^^), one usuahy considers the Poincare section (or stroboscopic 

map) ^{tq + nTy), with n = 1, 2, ... Because the Poincare section of the unperturbed system has a 
hyperbohc point at (1,0), the Poincare section of the perturbed system wih have a hyperbohc point 
of the same kind (saddle) in the vicinity of (1,0), provided e is small enough. An unstable (resp. 
stable) invariant manifold (resp.M^'^ ) will therefore exist in the vicinity of the hyperbolic point. 
If these two manifolds intersect transversally, they will have an infinity of such intersection points. 
The non-dissipative (area preserving) character of the 0{e^) system will induce huge stretching, 
and folding will inevitably follow. The phase portrait in the vicinity of the vortex will take the form 
sketched in figure EKb) which drastically differs from the non-chaotic case. Particle sedimentation 
will then be chaotic. Such intersection points can be detected by making use of the classical Melnikov 
method (see for example Guckenheimer & Holmes |8]) which consists in calculating the dot product 

where is the normal to the homoclinic trajectory at i^o(O) with (Yo,1^,~^z) 
right-handed, F (resp. E) is the intersection between this normal and (resp. ). If the two 
manifolds intersect transversally then (i(ro) will have simple zeros as tq varies. If (i(ro) remains 
strictly negative, then the relative position of the two manifolds will look like the one sketched in 
figure [2Kc) : lying outside and inside. In this case the particles located outside the cell will 
fall regularly without being catched into the cell, whereas the particles initially located inside the 
cell will spiral out and exit the cell. To order 0{e), the Melnikov "distance" d{To) is proportional 
to the Melnikov function[8] : 

/OO • / \ POO 

Voir) A \^}{^o{r)) sin ( -^(r + tq) ) dr - kV^ / Yo{t) A f o(r) dr 
-oo / J ~oo 

where tq is the starting time of the Poincare section of the dynamical system. By writing l^o = 
(^(T),?7(r)) and expanding the sine function, and noticing that ^ is an odd function, we are led to : 



M^i^o) = cos / i{T) sin -rp2 dr - kV^'j 

J — oo 



1=1 iiv- ilk)dT 



where 



is a purely numerical constant. As noticed above, both ^(r) and r/(r) are unknown, but they are 
purely numerical functions which can be determined from a numerical solution of equation ([8]). 
Finally, the Melnikov function of the homoclinic trajectory in the 1^ plane reads 

M^(ro) = cos (^ro) F - kVh (9) 

where F{s) is the sine transform of ^ and is a purely numerical function depending only on the shape 
of the initial vortex. Note that since ii) — r]^, is proportional to the curvature of the homoclinic 
trajectory, the constant part —kV^^ clearly manifests the contribution of a centrifugal effect due 
to particle inertia, as already observed for the onset of chaos in Stommel cells [1]. 

Because F{Q) = 0, the steady case a; = is straightforward : M^{tq) = —kV^^ < 0. This 
means that the manifold remains outside the cell, whereas remains in the inner side (like in 
figure [2](c)) : they will not intersect, and any particle released outside the cell will go down without 
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penetrating into the cell. Also, as already mentioned, particles released inside the cell will spiral 
out. We recover the fact that permanent suspension does not exist for such inertial particles in our 
flow, as already noticed by Wang & Maxey |12) . Rubin, Jones & Maxey |16) for other flows. 

In the unsteady case a; > 0, the Melnikov function has simple zeros if the amplitude of the 
oscillating term is larger than the constant term. The criterion for the appearance of chaotic 
particle sedimentation is therefore : 

F (^) > kVH (10) 

By making use of a numerical algorithm to solve Eq. ^ we obtain 7 ~ 45.8, together with F{s) 
which is plotted on figure [3l This function has a peak value max^ i^(s) ~ 1.14, so that no chaotic 
motion is expected to occur, under the present hypotheses, if (say) 

k-fV^ > 1.15 

In this case the Melnikov function remains strictly negative whatever the frequency of the pertur- 
bation : the particle motion is always regular. For kjV^ < 1.14 a homoclinic bifurcation can occur, 
provided w/V^ lies in an appropriate range ("chaotic window") like the one shown in figure [3j The 
particle motion is therefore highly affected by the oscillations of the vortex. 



3 Comparison with numerical solutions 

Figure H] shows the evolution of a particle cloud initially released inside the cell obtained by solving 
numerically the particle motion equation ([3]) with e = 0.1. The terminal velocity of the particles is 
Vt = 0.42 and the frequency of the perturbation is a; = 4V^ ~ 0.7. In addition k = 0.1, so that 
k'yV^ ~ 0.8 : according to figure [3] chaos is likely to occur in the vicinity of the separatrix. We 
indeed observe that the particle cloud is folded and stretched. Figure [5] shows two typical particle 
trajectories, initially released outside the cell for k'jV^ ~ 0.8. In the case 00 = 5V^ the motion is 
chaotic and the spatial length of the trajectory is larger than in the non-chaotic case (w = IIV^) 
because the particle is captured into the cell and spins there for a while. This property does not 
imply that the system is chaotic, but it is used in the following to detect the occurrence of chaos in 
a more systematic way, as proposed by Ziemniak &; Jung |20) . and as done in a previous paper [1]. 

Indeed, to check the predictive power of formula (fTO|) . we have run a set of computations where 
1000 particles are released slightly above the cell, with e = 0.2 and r = 0.02. The trajectory of 
each particle p is then calculated by introducing a random phase shift in the flow perturbation 
(simply replace cot by ut + Up in equation ([2|)). The calculation is stopped when the particle reaches 
a fixed bottom, below the cell. We then calculate the centred averaged path length : 

_ Ljuj) - L(oo) 
L[oo) 

where L{uj) is the average particle path. The quantity 6{uj) is plotted in figure [6] for four values of 
Vt- In addition, we have plotted the amplitude of the Melnikov function minus its constant part, 
that is F{ijj/V^) — k^V^, the positiveness of which implies chaotic particle motion. It appears that 
when F{uj/V^) - kjV^ < the average particle path length is closed to zero (i.e. all the particle 
paths have roughly the same length) , as expected if the two manifolds do not intersect (like in figure 
[2{c)). In contrast, a soon as F(uj/V^) — k^V^ > the average particle path length increases, and 
this manifests the fact that some particles have penetrated into the cell, as a consequence of the 
homoclinic bifurcation. 
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4 Discussion 



The calculations presented in this note show that a fixed vortex is sufficient to induce chaotic 
particle motion, under the sole effect of gravity and of the unsteadiness of the vortex. We have 
chosen to consider a fixed vortex with time-dependent intensity. One could argue that this choice is 
not realistic, since the intensity (circulation) of vortices is known to remain constant unless viscosity 
affects it (Kelvin's theorem). This is why we assumed that some appropriate boundary conditions, 
or an appropriate volumic force, was present to sustain the whole picture. Our goal being to show 
that gravity and unsteady differential rotation are sufficient ingredients to induce chaotic particle 
settling. Note also that the detailed shape of the vortex (here a rotation rate decaying like 
where r is the distance to the vortex centre) might not be of major importance, and that other 
decaying rotation rates could also lead to chaotic particle motion. 

In the present analysis, gravity plays a significant role since it is responsible for the appearance 
of a homoclinic trajectory in the leading-order dynamics. Like for particle settling in the vicinity of 
upward streamlines [Tj, particle inertia is opposed to the appearance of chaos, because of centrifugal 
effects, and tends to maintain the two invariant manifolds and away from each other. The 
flow unsteadiness, in contrast, tends to make these manifolds intersect. 

To check the predictive power of the Melnikov analysis we have computed the average length of 
particle paths, as proposed by Ziemniak & Jung |20) in the framework of fluid points trajectories in 
the wake of a cylindrical obstacle. These authors observed that the probability to flnd a trajectory 
whose path length increase is larger than some value s, is an exponential function of s. It would be 
therefore of interest to check whether such exponential distributions are visible also in the present 
case. Clearly, these statistics are linked to the area of the lobes (i.e. subsets of the phase space 
located between and W^, as shown in grey on flgure [2|^b)), which can be calculated from the 
Melnikov function. In addition, lobe dynamics can be used to calculate other global quantities 
like particle flux accross the separatrices (see for example Rom-Kedar, Leonard & Wiggins |15) . 
Balasuriya [5]). A detailed analysis of such integral quantities would therefore be of interest to 
characterize more precisely inertial particle transport in the chaotic regime. Further studies on this 
topic should be of interest. 

In the absence of particle inertia {vp = 0) the particle dynamics is always chaotic (just set k = 
in the Melnikov function ([9])) : this is a pure hamiltonian chaos, in the sense that the complete 
system is hamiltonian here, and very close to chaotic advection of perfect tracers. The Poincare 
sections in this case are very classical, and are shown in figure [7] in the case Vr = 1- As expected, 
KAM curves are visible outside the stochastic layer close to the homoclinic trajectory, indicating 
that some particles could be trapped for a while, and these curves are destroyed as the unsteadiness 
of the flow increases. 

The present mechanism could be applied to mixing devices in chemical engineering processes, 
where one could leave sedimentation act between every two stiring periods. To our knowledge, 
detailed mathematical analyses devoted to such devices have not been published so far. In another 
context, the stretch, sediment and fold mechanism could play a non- negligible role in the mixing 
of plankton (or the mixing of any other "particle") in the upper ocean. Indeed, under the com- 
bined effect of settling and of unsteady wind-induced (or temperature gradient induced) rotating 
flows, patches of sedimenting particles could perhaps be mixed efficiently. This point needs further 
discussions. 

The main conclusion of the present work is that an unsteady differential rotation is sufficient 
to induce chaotic heavy particle settling provided the still-fiuid terminal velocity of the inclusion 
is close enough to the flow velocity, and particle inertia is small enough. (In contrast, fluid point 
trajectories are very regular here, as fluid points go to and fro along portions of circle.) The 
mechanism, which can be called "stretch, sediment & fold", is sketched in figure [8] (for the sake of 
clarity, large amplitude oscillations are assumed there) : a particle cloud is stretched by a vortex. 
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then the vortex weakens and the cloud sediments, then the vortex starts again and the cloud is 
folded, and so on. This elementary mechanism is different from the one investigated by Vilela & 
Motter |18| with two blinking vortex-sources. Indeed, in the system investigated by these authors 
gravity is not a key ingredient for the appearance of chaos. The particle spirals out around the 
vortex-sources, until it reaches a limit cycle and remains suspended permanently. On this attractor 
the spiraling time of the inclusion is of the order of the period of the blinking. In our single- vortex 
case, it is the sedimentation time scale which has to match the blinking period. Moreover, these 
authors show that a cascade of period doublings occurs as the particle inertia decreases, leading to 
a strange attractor. In the present paper no attractor can be observed, since our system is reduced 
to a perturbed hamiltonian system. It could therefore be of interest to check whether, once this 
simplification is removed, a permanent suspension could appear in the vicinity of a single fixed 
singularity with time-periodic strength, and gravity. Further analyses should clarify this point. 
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Figure 1: Sketch of the streamhnes of a horizontal vortex (a) and of particle trajectories (to leading 
order e^) in the vicinity of this vortex (b). The homoclinic trajectory (dashed line) is attached to 
the saddle point A where the fluid velocity balances the terminal velocity. 
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Figure 2: Sketch of the Poincare section of the perturbed renormahzed system. In case (b) the 
manifolds and VF" intersect, leading to chaotic particle trajectories. In case (c) the two manifolds 
do not intersect : particles released inside the cell will spiral out, and those released outside will go 
round the cell. 
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chaotic sedimentation 

Figure 3: Plot of the amplitude of the Melnikov function, obtained by solving numerically the 
particle path over the homoclinic trajectory of the renormalized system. 
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Figure 6: Plot of the centred average path length {0{uj), black dots), together with the amplitude 
of the Melnikov function minus its constant part (solid line, formula (fTO|) ). As soon as this line 
is above zero, the separatrix gets broken, some particles therefore penetrate into the cell and 9{uj) 
increases. 
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Figure 7: Plot of 10 Poincare sections for = and Vr = 1- Because k = the criterion (jlOp 
obtained from Melnikov's analysis is always fulfilled, so that the stochastic layer in the vicinity 
of the separatrix always exists (except, of course, in the steady case £ = 0). The plots show the 
destruction of KAM tori as the amplitude of the perturbation increases. 
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(a) 1'^^ (c) (d) 




Figure 8: Sketch of the gravity-induced bhnking vortex effect, for large amplitude oscihations. 
Sketch (a) shows the initial particle cloud (black) and the vortex (white circle). In step (b) the 
particle cloud is stretched due to differential rotation. Then the vortex weakens and the cloud 
sediments (c). When the differential rotation restarts the cloud is folded. 
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